Amplitude squeezed fiber Bragg grating solitons 
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Quantum fluctuations of optical fiber Bragg grating solitons are investigated numerically by the 
back-propagation method. It is found for the first time that the bandgap effects of the grating act 
as a nonlinear filter and cause the soliton to be amplitude squeezed. The squeezing ratio saturates 
after a certain grating length and the fundamental Bragg soliton produces the optimal squeezing 
ratio. 



PACS numbers: 42.50.Lc, 42.65.HW, 42.65.Tg, 42.81. Dp, 05.45.Yv, 42.70.Qs 



m 

O 
O 

(N 
C 

m 

(N 



> 

o 
m 
o 

S3 

+-> 

G ■ 
cd 

3 : 
cr 



In the literature, various types of optical soliton 
phenomena have been studied extensively in the area 
of nonlinear optical physics. These include the non- 
linear Schrodinger solitons in dispersive optical fibers, 
spatial and vortex solitons in photorefractive materi- 
als/waveguides, and cavity solitons in resonators 0. 
It has also been well known that fiber Bragg grat- 
ings (FBGs) with Kerr nonlinearity can exhibit optical 
soliton-like phenomena known as the Bragg grating soli- 
tons 0, 0| ■ The FBGs are one-dimensional photonic 
bandgap crystals with weak index modulation. By uti- 
lizing the high dispersion of the FBGs near the band- 
edges, one can produce optical solitons in the anomalous 
dispersion side if the input pulse have suitable pulse- 
width and peak intensity. From the theoretical point of 
view, solitary waves in one-dimensional periodic struc- 
tures can travel with different group velocities and have 
been verified in some experiments [3j- Even for two- or 
three-dimensional nonlinear photonic bandgap crystals, 
solitary waves can also exist Q and have been observed 
recently 0- 

Most of the previous studies on Bragg grating solitons 
have been on the classical effects and there is almost no 
result on their quantum properties. The quantum theory 
of traveling- wave optical solitons has been intensively de- 
veloped during the past 15 years and several approaches 
have been successfully carried out to calculate the quan- 
tum properties of different traveling- wave optical solitons 
including the family of nonlinear Schrodinger solitons 
13, Q as well as the self-induced-transparency solitons 
8]. Bragg solitons belong to the class of bi-directional 
pulse propagation problems where the quantum theory 
is still lack of enough consideration. It is the aim of this 
study to bridge this gap by developing a general quan- 
tum theory for bi-directional pulse propagation problems 
and particularly applying the theory to the case of Bragg 
grating solitons. It will be shown that the output Bragg 
soliton pulses will quantum-mechanically get amplitude- 
squeezed and the squeezing ratio can be calculated theo- 
retically. 

In our modeling, we use the nonlinear coupled mode 
equations (NCMEs) to describe the two bi-directional 
waves propagating in a uniform FBG. We use the lin- 
earization approach to study the quantum effects of op- 



tical solitons in FBGs by extendingthe back-propagation 
method we previously developed [9j to the cases of non- 
linear bi-directional propagation problems. By following 
the same spirit of the back-propagation method, we will 
first derive a set of linear adjoint equations from the lin- 
earized NCMEs in such a way that any inner product 
between the solutions of the two equation sets is con- 
served during the time evolution. Under the lineariza- 
tion approximation, the measurements performed after 
the time evolution can also be expressed in terms of the 
inner product between the perturbed quantum field op- 
erator and a measurement characteristic function. By 
back-propagating the measurement characteristic func- 
tion to t = through the solution of the adjoint equa- 
tions, we can express the measured operator in terms 
of the input field operators which have known quantum 
characteristics. In this way, the variance of the measured 
operator as well as its squeezing ratio can be calculated 
readily for a given measurement characteristic function. 

To be more explicit, let us consider the wave propaga- 
tion problem in a one dimension fiber grating structure 
with the nonlinearity coming from the third order nonlin- 
earity of the optical fiber. With the self-phase modula- 
tion and cross-phase modulation effects, we model Bragg 
solitons by using the following NCMEs that describe the 
coupling between the forward and the backward propa- 
gating waves in a uniform FBG. 
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Here U a ,b{z, t) represent the forward and backward prop- 
agation pulses respectively, v g is the group velocity of the 
pulse, k is the coupling coefficient, As is the Bragg wave- 
length, 5 is the wavelength detuning parameter, and T 
represents the self-phase modulation coefficient. This set 
of NCMEs has approximate analytical soliton solutions 
for the case of infinite grating length, as is shown by 
Aceves and Wabnitz with the introduction of the mas- 
sive Thirring model Q- However, for gratings of finite 
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FIG. 1: Transmission (top) and optimal squeezing ratio (bot- 
tom) for Bragg solitons with different input intensities. 



length, no analytic solution can be found. So in our stud- 
ies we directly use the finite difference numerical simu- 
lation method with the parameters based on the first 



Vg dt\Ub) 

Vg dt\K J 

+ 



If we define the inner product operation according to 

</1I) = ~ / dz [f*J a + f a gt + r b 9b + fbfi] (5) 

and introduce the corresponding set of adjoint equations 
in Eq. (jlj, then it is easy to show that the inner product 
between the solutions of the two equation sets is pre- 
served along the time axis, 4z{ii \u) = 0. By taking ad- 



experimental report in the literature 3]. We consider 
a 60 ps FWHM sec/i-shaped pulse incidents into a uni- 
form grating with 15.0 cm -1 wavenumber detuning from 
the center of the bandgap. The coupling strength of the 
fiber grating is 10 cmT x . When the input intensity is 
below the required value for forming a solitary pulse in 
the FBGs (about 4.5 GW/cm 2 in this case), the peak in- 
tensity of the pulse will decrease along the propagation. 
On the other hand, when the input intensity is above 
4.5 GW/cm 2 , the peak intensity oscillates. Only when 
the nonlinearity can exactly compensate the dispersion 
induced by the FBGs, we can have a stable solitary pulse 
inside the grating. After obtaining these classical solu- 
tions, we now turn to the calculation of their quantum 
properties. 

Since for optical solitons the average photon number 
is usually very large, we can safely use the linearization 
approximation to study their quantum effects. By setting 

U a (z,t) = U a0 {z,t) +u a (z,t) 
U b (z,t) = Ubo(z,t) + u b (z,t) 

and substituting them into Eq. lllll'li for linearization, we 
obtain the linear quantum operator equations in Eq. @ 
that describe the evolution of the quantum fluctuations 
associated with the Bragg solitons. The quantum per- 
turbation fields u a (z, t) and u b (z,t) have to satisfy the 
following equal time commutation relations: 

[u a (z 1 ,t),ul(z 2 ,t)} = [u b (zi ) t),ul(z 2 ,t)] = 5(zi - z 2 ) 
[u a . b (z 1 ,t),u atb (z 2 ,t)} = [u\ b {zi,t),u\ h {z 2 ,t)] =0 



(3) 



(4) 
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vantage of this property, we can express the inner product 
of the output quantum perturbation operator with a pro- 
jection function in terms of the input quantum operator. 
This will allow us to calculate the quantum uncertainty 
for the inner product of the output quantum operator 
with any given projection function. For the photon num- 
ber measurement, the pulse shape of the projection func- 
tion is simply the classical output solution 10]and the 
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FIG. 2: Squeezing ratio for Bragg solitons with different 
length of FBGs. 



squeezing ratio can be calculated according to: 

£/ T s = var[{f(z)\u(z,t = T)} 
var[(f(z)\u(z, t = 0)] 
var[(F T (z)\u(z,t = 0)] 
var[{f(z)\u(z,t = 0)] 

Here f(z) is the original projection function and Ft(z) 
is the back-propagated projection function. 

In the following we shall present our results of simula- 
tion. The transmission of the FBGs with different input 
intensities for a constant FBG length (50 cm) is shown 
in the top curve of Fig. ^ The calculated photon num- 



ber squeezing ratio is shown in the bottom for the same 
parameters. The squeezing ratio decreases monotonically 
when the intensity is smaller than the intensity of the fun- 
damental soliton but will begin to oscillate strongly when 
the intensity is much larger than that of the fundamental 
soliton. The oscillation behaviors of the FBG transmis- 
sion and the squeezing ratio match very well. That is, 
the squeezing ratio has a local minimum when the trans- 
mission has a local maximum. Intuitively the periodic 
grating structure acts like a spectral filter which can filter 
out the noisier high frequency components in the soliton 
spectrum and produce a net amplitude squeezing effect 
just as in the previous solution amplitude squeezing ex- 
periments where a spectral filter is cascaded after a non- 
linear fiber 0, • It is also intuitively clear that larger 
amplitude squeezing should occur when the transmission 
curve is saturated. Fig. [3] shows the dependence of the 
calculated squeezing ratios for different FBG lengths. If 
we only consider the gratings with the length longer than 
1 cm, we find that the squeezing ratio monotonically de- 
creases with the FBG length and saturates at the length 
around 60 cm. Intuitively this is because the filtering ef- 
fect of the grating will unavoidably introduce additional 
noises on the light fields and eventually cause the squeez- 
ing ratio to become saturated. 

To summarize, we have developed a general quantum 
theory for bi-directional nonlinear optical pulse propaga- 
tion problems and have used it to calculate the squeez- 
ing ratio of Bragg solitons in one-dimensional photonic 
bandgap crystals for the first time. We find that the out- 
put pulses can get amplitude squeezed and the squeezing 
ratio exhibits interesting relations with the fiber grating 
length as well as with the intensity of the input pulse. It 
will be very interesting to see if one can actually observe 
these effects experimentally in the future. 



[1] For a review of various soliton phenomena, please see 
"Spatial Solitons", edited by S. Trillo and W. Torruellas 
(Springer, 2001). 

[2] A. B. Aceves, and S. Wabnitz, Phys. Lett. A 141, 37 
(1989). 

[3] B. J. Eggleton, R. E. Slusher, C. M. de Sterke, P. A. 

Krug, and J. E. Sipe, Phys. Rev. Lett. 76, 1627 (1996). 
[4] A. A. Sukhorukov and Y. S. Kivshar Phys. Rev. Lett. 

87, 083901 (2001). 
[5] J. W. Fleischer, M. Segev, N. K. Efremidis, and D. N. 

Christodoulides, Nature (London) 422, 147 (2003). 
[6] P. D. Drummond and S. J. Carter, J. Opt. Soc. Am. B 



4, 1565 (1987). 
[7] Y. Lai and H. A. Haus, Phys. Rev. A 40, 844 (1989); 

Phys. Rev. A 40, 854 (1989). 
[8] Y. Lai and H. A. Haus, Phys. Rev. A 42, 2925 (1990). 
[9] Y. Lai and S.-S. Yu, Phys. Rev. A 51, 817 (1995). 
[10] H. A. Haus and Y. Lai J. Opt. Soc. Am. B 7, 386 (1990). 
[11] S. R. Friberg, S. Machida, M. J. Werner, A. Levanon, 

and T. Mukai, Phys. Rev. Lett. 77, 3775 (1996). 
[12] T. Opatrny, N. Korolkova, and G. Leuchs, Phys. Rev. A 
66, 053813 (2002). 



